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ABSTRACT 

Long-term evolution of a stellar orbit captured by a massive galactic center via suc- 
cessive interactions with an accretion disc has been examined. An analytical solution 
describing evolution of the stellar orbital parameters during the initial stage of the 
capture was found. Our results are applicable to thin Keplerian discs with an arbitrary 
radial distribution of density and rather general prescription for the star-disc interac- 
tion. Temporal evolution is given in the form of quadrature which can be carried out 
numerically. 
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1 INTRODUCTION 

Evolution of the orbit of a star under the influence of inter- 
actions with an accretion disc has been studied by numerous 
authors because this situation is relevant to inner regions of 
active galactic nuclei. The trajectory of an individual star is 
determined mainly by gravity of the central mass and sur- 
rounding stars while periodic transitions through the disc 
act as a tiny perturbation. The final goal is to understand 
the fate of a star, transfer of mass and angular momentum 
between the star and the disc, and also to determine how 
star-disc interactions influence the distribution of stellar or- 
bits near a massive central object. An important and difficult 
task is to estimate the probability that a star gets captured 
from an originally unbound orbit, and to determine whether 
this probability is significant compared to other mechanisms 
of capture. 

Orbital evolution of a body crossing an accretion disc 
has been discussed with various approaches, first within the 
framework of Newtonian gravity, both in theory of the solar 
system (Pollack, Burns & Tauber 1979; McKinnon & Leith 
1995) and for active galactic nuclei (Goldreich & Tremaine 
1980; Syer, Clarke & Rees 1991; Artymowicz 1994; Podsi- 
adlowski & Rees 1994). These studies have been general- 
i zed in order to account for the effects of general relativity 
( Vokrouhlicky & Karas 1993 ) and to model a dense star clus- 
ter in a galactic nucleus (Pineault & Landry 1994; Rauch 



1995). It has been recognized that detailed physical descrip- 



tion of the star-disc interaction is a difficult task (Zurek, 
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Siemiginowska & Colgate 1994, 1996). In this letter a sim- 
plified analytical treatment of stellar orbital parameters is 
presented, describing the initial stage of star-disc collisions 
(when the star crosses the disc once per revolution). A great 
deal of our calculation is independent of microphysics of 
star-disc interaction. It is shown how our solution matches 



the corresponding Rauch's (1995) solution which is valid in 
later stages, when eccentricity of the orbit becomes small 
enough. 

In the next section our approach to the problem is 
formulated and an analytical solution is given. Then, in 
Sec. 2.2, further details about the derivation of the results 



are presented, and finally a s imple example of the orbital 
evolution is shown in Sec. 



2.3 



2 STELLAR CAPTURE BY A DISC 
2.1 Formulation and results 

Newtonian gravitational law is assumed throughout this pa- 
per. Our solution is based on the following assumptions: 

(i) the disc is geometrically thin and its rotation is Keple- 
rian; 

(ii) at the event of crossing the plane of the disc, velocity 
of the star is changed by a tiny quantity. This impulse is 
colinear with the relative velocity of the star with respect to 
the material forming the disc; 

(iii) the star crosses the disc once per revolution (the model 
is applicable to the initial phase of the stellar capture). 

The first item is a standard simplification in which the disc is 
treated in terms of vertically integrated quantities, while the 
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second one can be expressed by the formula for an impulsive 
change of the star's velocity: 



Sv = E(r,v)v le i: 



(1) 



E is an unconstrained function, given by a detailed model of 
the star-disc interaction, and v re \ is relative velocity of the 
star and the disc material. We stress, that our results are 
uniquely based on this assumption of colinearity, Sv oc v le v, 
the coupling factor E is arbitrary and it can be as complex as 
necessary. In particular, E contains information about the 
surface density i of the disc (i = outside an outer edge 
r = Rd of the disc). The form of E must be specified only 
for examination of temporal evolution of orbital parameters. 
We will assume, in analogy with Rauch (1995), a simplified 
formula for 



E(r,v) 



* (r) — £ . 

v± 



i + 



In A, 



(2) 
(3) 



when it is needed for purpose of an example. In eq. (|2j) /?* 
denotes radius of the star, m* is its mass, escape velocity 
(v 2 = 2Gm*//?*) ; v± is normal component of the star's 
velocity to the disc plane, and In A is a usual long-range 
interaction factor (Coulomb logarithm). 

The star's orbit is traditionally characterized by the Ke- 
plerian osculating elements: semimajor axis a, eccentricity 
e, inclination I to the accretion disc plane, and longitude of 
pericenter ui (measured from the ascending node). A derived 
set of parameters turns out to be better suited to our prob- 
l/a, n_ 



lem: 



/'• 



cos/, and k = ecostj. We 



will show (Sec. 2.2) that evolution of a stellar orbit following 
the capture by a disc can be written in a parametrical form: 

«(C) = ^(0{«o^- 1 (Co) + a 2 [V(C)-^(Co)]} , (4) 
77 2 (C) = C<MC) W-%)a- 2 + [<MC)-^(Co)]} , (5) 

tic) = Vc + o(o, (6) 

l*(C)l = C-i, (7) 

where the auxiliary functions 0(C), 0(C) an d ip(C) read 
^ ( 1 + * h - r. -x- r.t* 
l 



<MC) = 
0(0 = 
iKO = 



c 



C 



l ± vi - c + cc 



2 + 



i - c + cc 
c 



i ± vi - c + cc 



(8) 
(9) 
.(10) 



Formal parameter C of the solution decreases from its initial 
value Co = 1 + eo | coswo | to the final value £f, given by 



Cf 



2R d a 2 
1 + i? d a 2 



(11) 



At this instant, the orbit starts crossing the disc twice per 
revolution and our solution ceases to be applicable. Obvi- 
ously, integration constants in (§)-(0) are determined in 
terms of the initial Keplerian orbital elements (ao, eo, Io, wo) 
by 



ao = 
Co = 
a = 

C = 
where 



fflo 

1 + eo | cos ujq | 

Co ' 
zo(zo 



2) 



(^o + l) 2 1-Co 



zo = — - 



cos/* = 



Co 



(17) 



(12) 
(13) 
(14) 

(15) 
(16) 

VCo(cos/ - VCo) 

Upper signs in (^)-(|Io|) are for the initial inclination Io 
greater than a critical value /* given by 

1 

TTo' 

lower signs apply otherwise. Integration constant C is singu- 
lar (C — > oo) for Io = /* (zo = — 1), and the solution can be 
simplified further. For instance, /i(C) = l/\/C f° r an values 
of C down to Cf ■ 

Solution (^)-(Q) can be extended easily to the case of 
initially parabolic orbits by setting a?o = 0, eo = 1 and 
a 2 — (Co/2-Rp). Here, R p denotes pericenter distance of the 
initial parabolic orbit. 

It is worth mentioning that the parameter C does not 
determine the time-scale on which the evolution takes place. 
Indeed, eqs. (^)-(Q) do not provide temporal information be- 
cause it depends on the precise form of E in eq. (Q). On the 
other hand, the strength and the beauty of the parametric 
solution (Q)-(Q) is in its independence on a particular model 
for E. We will also illustrate an example of temporal evolu- 
tion later in the text, and only for this purpose the form of 
E will be needed. Assuming relation (^)-Q one obtains 



t-to = 



ct 3 VGM E c 



Co 



(18) 



where to is initial time, M is the central mass, and E c = 
(7r/? 2 /m*) i (r c )C with r c = a~ 2 being radial distance of the 
point of intersection with the disc. Function v = v-cexjv^ is 
determined by orbital parameters which themselves depend 
on C according to eqs. (|4])-(^) 



We note that Rauch (1995) conjectured that function 
R = arf cos 4 (//2) remains nearly conserved along the evolv- 
ing stellar orbit and he used this function for estimates of 
the radius of the final, circularized orbit in the disc plane. 
In our notation, 

R{0 = 4^ [ 1 + V / C + ^(0] 2 - (19) 

Hereinafter, we show that R(C) is a well-conserved quan- 
tity at later stages of the orbital evolution (when eccentric- 
ity is sufficiently small), but it fails to be conserved at the 
very beginning of the capture when the orbit is still nearly 
parabolic, close to an unbound trajectory. 



2.2 Details of the solution 

In this section, we present more details of the derivation of 
the solution given above. 

Taking into account the fundamental assumption (Q), 
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Figure 1. Inclination 7(f) of a captured stellar orbit (measured 
in units of initial inclination Jo) vs. Jo. This graph corresponds to 
initially parabolic orbits with pericenter distance equal to the disc 
radius. The curves are parametrized by f . Temporal evolution of 
some particular orbit starts with I = Io, f = 2, and goes down 
along the vertical line, across f = const curves. Our analytic 
solution (solid lines) is valid in the region 2 > £ > 1, and it 
corresponds to eccentric orbits. Circularization time-scale is equal 
to the time-scale necessary for grinding the orbit to the disc plane 
if Jo < J* = 45° while the former is shorter than the latter for 
orbits with Jo > J*. At £ = 1, i.e. zero eccentricity, our solution 
matches Rauch's R = const solution (dashed lines). 



one can easily demonstrate that the Keplerian orbital ele- 
ments are perturbed at each transition (due to interaction 
with the disc material) by quantities 



S(V^) = 

5{yfari) = 

5{-s/ar)[i) = 

5(k) = 



E ^fa-q fi , 
E Var] /2 , 
E Varj fa , 
2E(l + *)/ 2 . 



Here, we introduced auxiliary functions 



/'• 



h = V 

h = i- 

h = v 7== ■ 

Combining eqs. ( ^l[ ) and 



{(1 +k)[l- vVTTk] + e 2 + k) 



1 



(20) 
(21) 
(22) 
(23) 

(24) 
(25) 

(26) 



i) we find that 



(27) 



VT+k = 

is conserved at the star-disc interaction. Hence, a is con- 
stant whatever the evolution of elements a, e and k is. In 
fact, condition ( gjj ) states that the initial Keplerian orbit has 
the same radius of intersection as the final orbit, after suc- 
cessive interactions with the disc. Longitude of the node is 
also conserved and can be set to zero without loss of general- 
ity. The above-given formulae (pii|)-(p7[) correspond to k > 
(i.e. \uj\ < it/2); for k < one should replace k — > —k. We 
note that all these relations can be easily reparametrized in 
terms of binding energy E = GM/ (2a), angular momentum 
L = V GMarj, and component of the angular momentum 
with respect to axis, L z = V GMar/fj,. 



Combining eqs. @ and @ with the help of Q, and 
introducing auxiliary variables y = \/arjfi and x = ^/arj, we 
obtain differential equation 



dy _ x(ay — 1 

dx (75 



(28) 



(We were allowed to change variations, 8, to the differentials, 
d, assuming an infinitesimally small perturbation of the stel- 
lar orbit at each intersection with the disc.) The Abel-type 
differential equation ( |28| ) can be solved beautifully by stan- 
dard methods of mathematical analysis (see, e.g., Kamke 
1959). An appropriate change of variables gives directly a 
solution for the evolution of inclination, eq. ([]). 

Similarly, considering ( |2c| ) and ( ^l| ) in terms of a = 1/a, 
we obtain, after brief algebraic transformations, 



\/C0(C)^ + «(C) 



C + y/ce (0 



(29) 



with #(£) given by eq. (g). This is a linear differential equa- 
tion, integration of which yields «(C) an d then, by eqs. (|J)- 
(§), also 77(C). 

At this point, one can see that Rauch's (1995) "quasi- 
integral" R is changed at each passage across the disc ac- 
cording to 



(30) 



S(lnR) = 2E 1 



Realizing that k ~ e we conclude that In J? indeed stays 
nearly constant at later stages of the orbit evolution, when 
eccentricity has decreased enough. On the other hand, at the 
very beginning of the capture process, when eccentricity is 
still high, R fails to serve as a quasi-i nteg ral of the problem. 
Instead, its evolution is given by eq. QlSh. 

For temporal evolution, eqs. (J2C|)— (p7|) must be supple- 
mented by additional relation, 

5(t) = -^a 3/2 , (31) 



which determines interval between successive intersections 
with the disc. Combining eq. ([H]) with (|2?]) one obtains 
separated differential equation which yields formula (|l8[). 
Recall that this last step requires assumption (Q) about the 
form of E. In the present case, 



2C(i - C - VCO) + C - v 2 



(32) 



1 - C - 2^(0 - 2 

Relation for time is apparently too complicated to be inte- 
grated analytically but numerical evaluation is straightfor- 
ward. 



2.3 Example 

We shall briefly demonstrate some properties of the analyt- 
ical solution from Sec. 2.1. 

We examine parabolic orbits (ao = 0, eo = 1) with peri- 
center in the disc plane (uo = 0) , and the pericenter distance 
R p equal to the disc radius (R p — iid = r c ). Initial inclina- 
tion Io of the orbit to the disc plane is a free parameter in 
this example. Evolution of this set of orbits is split into two 
phases. 

First, we let the orbits evolve according to the solution 
of eqs. (M-(pj) from the initial value £o = 2 of the formal 
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Figure 2. Function R (normalized to its initial value Ro) vs. 
initial inclination io of the orbit for different values of parameter 
f. The ^-dependence is stronger for orbits with small To- Once the 
orbit becomes circular with f = 1, R reaches its terminal value 
and does not change any more. Therefore, R does not acquire 
values inside the shaded region. 



Figure 3. Time i c (C) °f orbital circularization of parabolic orbits 
with initial inclination Io, as in Fig. 1. Here, time (arbitrary units 
on vertical axis) is recorded starting from eccentricity e = f — 1 
(given with each curve) down to e = 10 — 3 (nearly circular orbit). 
Notice the change in form of the curves at critical inclination 
h = 45? 



parameter £ to its final value Q = 1. Figure |l| illustrates 
the evolution of the inclination 7(C), measured in terms of 
the initial value Jo. Notice that the critical inclination of 
eq. ([n]) is 45? Eccentricity of the orbits under consideration 
decreases according to a simple formula e(£) = £ — 1 (inde- 
pendently of Io), leading eventually to circularized orbits at 
C = Cf ■ We observe that orbits with Io < I* terminate at the 
final state If = 0, suggesting that the circularization time- 
scale is comparable to that necessary for grinding the orbit 
into the disc plane. On the other hand, when Jo > I* the 
final circular orbits remain inclined significantly to the disc 
plane. (/ > 90° corresponds to retrograde orbits.) Hence, 
for those orbits the circularization time-scale is shorter than 
the time necessary for tilting the orbit to the disc plane. 
Additional time to incline a circularized orbit is not much 
longer than circularization time, however The difference is 
typically a factor of 10 for highly retrograde orbits. 

Secondly, we examine the evolution of circularized or- 
bits which started with Jo > /* and have settled to nonzero 

= 1). Because these orbits have zero eccentricity, there 
exists Rauch's integral in the form R\ = a 2 acos 4 (I /2) = 
£ cos 4 (1/2). Here, we adopted a formal continuation of the 
£ parameter to values smaller than unity (in this phase, 
£ = a/R c \). For each orbit, we calculate the value of Ri = 
R(<^ = 1), so that the inclination is given by 



(33) 



Obviously, a given orbit terminates its evolution at £ = 4R\ 
when it is pushed completely into the disc plane. Dashed 
curves in Figure |l| correspond to constant values of £ < 1. 

Figure □ illustrates how function R(C) changes during 
the first circularization phase of the evolution. For each orbit 
we have chosen the same steps in £ (0.2) in the range 1.8 > 

t We thank the ref eree for pointing out this fact, confirmed a lso 
by other estimates ( [gycr ct al. 1991 ; McKinnon fc Lcith 1995| ) . 



£ > 1.0, as in Figure^, and we computed the corresponding 
values of R(£) from eq. ( |l9| ) . Our results agree with Rauch's 
( 1995| ) finding, namely that R(£) is conserved up to a factor 
of w 2 for orbits with large eccentricity. During the second 
phase of the evolution the i?-function is constant. 

Figure ^ shows time intervals i e (C) which elapse in the 
course of gradual circularization when eccentricity decreases 
from e = C — 1 to some terminal value (here, terminal ec- 
centricity has been fixed to e = 10 -3 ; notice that t e goes 
to infinity for terminal eccentricity e — » 0). We have ver- 
ified the graph also by direct numerical integration of the 
corresponding orbits. Numerical factor standing in front of 
integral on the right-hand side of eq. (^) can be written in 
physical units in the form 



10' 



10 3 R e 



9/4 



Re. 



Go • 



10 3 Ji g 
R* 



^ yrs; (34) 



R e = 2GM/c? and R gir = 2Gm±/c 2 axe gravitational radii of 
the central mass and the star, respectively. A typical surface 
density profile of the disc has been assumed, as in eq. (1) 
of Rauch (1995). The value of £ ~ 10 3 corresponds to the 
estimate in addendum to Zurek, Siemiginowska & Colgate 
(1996). 



3 CONCLUSION 

We have found a solution describing the evolution of orbital 
parameters of a star orbiting around a massive central body 
in a galactic nucleus and interacting with a thin Keplerian 
disc. The solution is in a parametrical form valid for an ar- 
bitrary radial distribution of density and a very broad range 
of models of the star-disc interaction. Temporal evolution 
can be given in terms of quadrature provided the star-disc 
interaction is specified completely (in terms of function E). 
Our approach can be applied to other situations but the 
form of eq. (Eq) is linked to the assumption about interac- 
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tions, eq. ([j]). Also the situation when the orbit intersects 
the disc twice-per-revolution requires a specific form of E 
to be given and, most likely, it does not allow a complete 
analytic solution. 

Our solution thus describes the initial phases of the stel- 
lar capture (large eccentricity) and it matches smoothly the 
low-eccentricity approximation. Apart from an interesting 
form of analytical expressions, our approach is useful as a 
part of more elaborate calculations. In an accompanying de- 
tailed paper, additional effects are taken into account (e.g., 
gravity of the disc) and distribution of a large number of 
stars is investigated (Vokrouhlicky & Karas 1997, submit- 
ted to MNRAS). 

We thank the referee for comments concerning tempo- 
ral evolution of the orbits and for other suggestions which 
helped us to improve our contribution. We acknowledge sup- 
port from the grants GACR 205/97/1165 and GACR 202/ 
96/0206 in the Czech Republic. 



REFERENCES 

Artymowicz P., 1994, ApJ, 423, 581 

Goldreich P., Tremaine S., 1980, ApJ, 241, 425 

Kamke E., 1959, DifTerentialgleichungen: Losungsmethoden und 

Losungen (Leipzig) 
McKinnon W. B., Leith A. C, 1995, Icarus, 118, 392 
Pineault S., Landry S., 1994, MNRAS, 267, 557 
Podsiadlowski P., Rees M. J., 1994, in: The Evolution of X-ray 

binaries, eds. S. S. Holt & C. S. Day, (AIP Press, New York), 

p. 71 

Pollack J. B., Burns J. A., Tauber M. E., 1979, Icarus, 37, 587 
Rauch K. P., 1995, MNRAS, 275, 628 

Syer D., Clarke C. J., Rees M. J., 1991, MNRAS, 250, 505 
Vokrouhlicky D., Karas V., 1993, MNRAS, 265, 365 
Zurek W. H., Siemiginowska A., Colgate S. A., 1994, ApJ, 434, 46; 
1996, ibid, 470, 652 



© 0000 RAS, MNRAS 000, 0-| 



